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Abstract 

We obtain exactly the vacuum expectation values ( (dtp) 2 (dtp) e lOL(p ) in the sine-Gordon 
model and ( L_ 2 L-2&Lk ) in $1,3 perturbed minimal CFT. We discuss applications of these 
results to short-distance expansions of two-point correlation functions. 
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1. Introduction 



One-point Vacuum Expectation Values (VEV) of local fields are important character- 
istics of Quantum Field Theory (QFT) vacuum. Operator Product Expansions (OPE) [|l| 
give rise to short-distance expansions for multipoint correlation functions which involve 
the one-point VEV as the basic ingredients 0,0]. At the same time the one-point VEV 
are nonperturbative objects and no systematic techniques for their evaluation is known. 
Some results for these quantities are available from numerical analyses (see e.g. || for 
numerical results in 2D QFT). Recently some progress has been made in evaluation of the 
one-point VEV in 2D integrable QFT |||6|,0]. In these papers the one-point VEV of the 
primary fields in some integrable QFT, including the sine-Gordon model and $i ; 3 per- 
turbed minimal CFT, were found exactly. On the other hand complete characterization of 
the correlation functions requires the knowledge of the VEV of all local fields, including 
the descendant operators. In the present paper we address the problem of calculating the 
one-point VEV of descendant fields. It was shown in [[|[7| that the VEV of the primary 
fields in sine-Gordon model (and in similar integrable QFT) satisfy remarkable "reflection 
relation" which involves the "reflection S- matrix" of Liouville CFT || , and their one-point 
VEV can be obtained as appropriate solutions to these relations. We will show here how 
this approach can be extended to the descendant fields and explicitly evaluate the VEV 
of the simplest nontrivial descendants in the sine-Gordon model and in $13 perturbed 
minimal models. 

Let us introduce basic notations and state main results of this work. The sine-Gordon 
model is defined by the Euclidean action 

A SG = j d 2 x I y^(cM 2 - 2fx cos(^) j , (1.1) 

where \i and (3 are parameters, < /? < 1. The simplest local fields in this QFT are 
the exponentials e lOLlf . Exact VEV ( e lOL{p )sg of these fields are found in ||. Here we will 
consider more general local fields of the form 

(d ni (p) (<9 n V) . . . (<9 nj » (<9 m V) (<9 m V) . . . {B mK <p) e iaip , I fte a I < J- , (1.2) 

where d = d z , d = d s and z,z are complex coordinates, z = X\ + ixi, z = X\ — 1x2- 
Precise definition of these fields in ( |1 . 1|) requires specification of their renormalizations. 
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We adopt the scheme in which the renormalized fields (|1.2| ) have definite scale dimensions 
(see e.g. []9|||]). Some details are presented in Sect. 2. In Sect. 3 we generalize the "reflection 
relations" of |6|,[7] to the fields ( |L2"| ) and use these relations to obtain the one-point VEV 
of the simplest nontrivial field of this kind, 



( (d^ (dip) ) SG = -( ) SG 



AI 



r(!) 



X 



r(-f + !)r(f + |)r( 



2 /3 2^ M 



I _u 21 

2^/3 



§) 



(1.3) 



where 



r(i + f-I)r(i-f -|)r(l + f + |)r(f-f + |) 

< " 2 



1-/3 2 

and M is the sine- Gordon soliton mass which relates to the parameter \x in as NTT 



1.4) 



res 5 



7rr(l-/? 2 ) 



M 



V^r(i + f) 
2 r(|) 



£\ -i 2-2/3 2 



(1.5) 



Let us quote here the simpler form the expression ( |1.3| ) assumes in the case a = El, 

( (dip) 2 (dip) 2 ) SG = -7T 2 M 4 tan 2 «/2) . (1.6) 



The sine-Gordon QFT is closely related to the Minimal CFT [|12] perturbed by the 
operator $i 3 , i.e. the QFT defined by the action 



M p/p , + X J d 2 x $i, 3 (x) 



;i-7) 



where Ai p / P ' stands for the formal action of the minimal model. Here we consider the 
massive case A > 0. As is well known, this CFT can be obtained from ( p..l|) with 
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p> — p 



1 It is interesting to note that Eq.(1.6) implies remarkably simple relation between the VEV of 
different fields associated with the sine-Gordon energy-momentum tensor T^ v . Denoting as usual 
T = T zz , T = Tzz and G = T Z z = jT" the irreducible spin components of and using the 



known result for (Q)sg [H| we have 
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by quantum group restriction 1 1^,14]. This relation was used in |J|j]] to obtain the VEV 

in (Q) 



of all primary fields $2^ in the QFT (|1 . 7|) . In a similar fashion, the fields (|1.2|) 
are related to the descendant fields in (|T77|) . In particular, the result (|1.3| ) is sufficient to 
derive the expectation values of the descendant fields L_2-L_2$;,A;, namely 
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Here | S ), s = 1,2, . . . ,p — 1 is any one of p — 1 degenerate ground states of the QFT 
( |1 . 7|) , and M is the mass of its fundamental kink. This mass is related to A in (|1.7|) as |L0] 
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(1.9) 



In fact, the "reflection relations" admit certain ambiguity to their solution, akin to the 
"CDD ambiguity" in the factorizable S-matrix theory (see e.g. |[L5|1 ). To some extent the 
ambiguity can be narrowed by taking into account the "resonance conditions" (see Sect. 2). 
In (|1.3| ) we have fixed this ambiguity by choosing the "minimal solution" - the simplest 
solution compatible with the resonance conditions (choosing the minimal solution is a 
common practice in the S-matrix theory). This choice is confirmed in Sect. 4-5, where ( |1.3j ) 
is checked against results obtained in (|1.1|) using semiclassical approximation (Sect. 4) and 
ordinary Feynmann perturbation theory (Sect. 5). Moreover, the special case ( |1.6| ) can be 
obtained directly from exact lattice theory of the XYZ model, as we show in Sect. 6. Finally, 
in Sect. 7 we use the Eq.( J1.8[ ) (more precisely, its particular case I = k = 0) to extend by 
one more order the short-distance expansion of the two-point correlation function of the 
Scaling Lee- Yang Model @. 



2. Descendant fields and Operator Product Expansions 



The sine-Gordon model ( |1.1| ) can be regarded as Gaussian CFT 



•AGai 



d 2 x (d u Lp) 2 
Ibn 



(2.1) 



perturbed by the relevant operator 2 cos(/3<^). Let Tcauss be the space of local fields in the 
CFT ( |2.1| ). This space is spanned by the fields (|1.2j). In the free field theory ( |2.1| ) these 



composite fields are defined through usual Wick ordering with 

( <p(z, Z) ip(0, 0) ) Ga uss = -2 log(zz) . 

With this definition the field ( |1.2| ) has the dimensions (A, A) = (a 2 + 1, a 2 + 1 ), where the 



integer "levels" /, / are the sums / = J2iLi n ii I = ^2f=i m j- The 



3- 



sum 



D = A + A = 2a 2 + l + l (2.2) 



is the scale dimension of the field (|1.2|) while the difference S = I — I coincides with its 
spin. Note that some linear combinations of the fields ( |1.2|) are total derivatives of other 
local fields, for example 

-a 2 {d<p){5<p) e ia<p = dde iaip , (2.3) 
{d 2 ip){dip) 2 e iaip + ia(d(p) 2 (3(p) 2 e* a(f = d((dip)(dip) 2 e iav ) , 
(d 2 ip)(8 2 ip) e ia ^ + ia(dip) 2 (d 2 ip) e ia ^ = d((d<p)(5 2 <p) e iaip ) . 
These elementary relations follow from the equation of motion of ( |2.1| ), 



(2.4) 



dd^p = . (2.5) 

The perturbation in (|1 . 1| ) leads to additional divergences in the matrix elements of the 
fields (|1.2[) , and so certain counterterms have to be added to ( |1.2| ) in order to compensate 
for these divergences. These counterterms contain local fields of the same spin, with cutoff- 
dependent coefficients. However, as long as the perturbation is relevant (i.e. < (3 2 < 1), 
the situation is relatively simple. For a given field (|1.2j ) it suffices to add only the fields 
of lower scale dimensions. In particular for a given field (|1.2|) there are finitely many 
counterterms. As usual these counterterms are not completely determined by the sole 
requirement that they absorb all the divergences, there is always a possibility to add 
finite counterterms. However, except for when certain resonance conditions are satisfied 
(see below), the ambiguity is eliminated completely if one demands that the resulting 
renormalized field (i.e. the unperturbed field ( |1.2[ ) plus the counterterms) has definite 
scale dimension, which then coincides with ( |2.2| ). We say that the field Oi has n-th order 
resonance with the field Oj if the dimensions of these fields satisfy the equation 



Di = Dj + 2n(l- (3 2 ) 
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(2.6) 



with some positive integer n. If this resonance condition is satisfied an obvious ambiguity 



Oi -> O l + const n n Oj (2.7) 

in defining the renormalized field Oi with the scale dimension Di typically results in the 
logarithmic scaling of Oi. 

We define the field ( |1.2|) in the perturbed theory (|1 . 1|) as the renormalized field sat- 
isfying the following conditions: i) ( |1.2| ) has definite scale dimension ( |2.2| ) and ii) the 
short-distance limit of its correlations coincides with the correlations of the corresponding 
field ( |1.2j ) in CFT ( |2.1|) . As explained above, in non-resonant cases this definition is unam- 
biguous. In resonant cases it is not complete and one has to impose additional conditions 
to fix the ambiguity (|2.7|) . For a given field ( |1.2| ) the resonances ( |2.6|) occur at isolated 



values of a. In this paper we are interested in generic values of this parameter and so we 
will not elaborate any specific convention concerning resonant cases. It suffices to note 
that if the matrix elements of (|1.2| ) are viewed as the functions of a the resonances show 
up as the poles in this variable. 

It is important to note that under this definition of the fields (|1.2|) in ( |1 . 1| ) these fields 
satisfy exactly the same relations ( |2.3| ), (|2.4| ) (as well as similar higher-level relations) as 
the fields in (|2.1| ). In this sense the symbol ip in ( |1.2|) is the subject to the "free" equation 
of motion (|2.5| ) rather then to the full equation of motion of the sine-Gordon theory. This 
is not a contradiction because it is rather the unrenormalized fields in (|1 . 1| ) that satisfy 
the full equation of motion; the renormalized ones differ from those by counterterms. In 
this paper we are interested in the one-point VEV of the fields (|1 . 2|) in the QFT (|1 . 1[) . Let 
us mention here some elementary properties of these VEV. First, only the fields of zero 
spin, i.e. with / = /, have nonvanishing one-point VEV. Next, the one-point VEV of the 
fields which are total derivatives vanish. According to (|2.3| ) , ( |2.4| ) the only spinless field of 
the level / = 1 has vanishing VEV, 

ad<p)(d<p)e io "<>)g G = 0, (2.8) 

There are four linearly independent spinless fields on the level 1 = 2. However the following 
relations among their VEV are simple consequences of 



(0V) {d 2 ip)e^ ) SG = -ia ( {dip)\d 2 ip) e ia * ) SG = 

-ia((d 2 ^(dcp) 2 e^)sG = -a 2 ((d<p) 2 (d<p) 2 )sg • ^ 
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These relation allow one to express all / = 2 VEV through the VEV 



{{dip) 2 {dip) 2 e^) SG . (2.10) 

Starting from the level / = 3 there are additional "kinematic" relation among the VEV 
following from the existence of higher local Integrals of Motion in the QFT (|1 . 1|) |fl6f , but 
we will not discuss them here. Instead we will concentrate attention on the VEV ( |2.10|) . 



Let us make here a simple remark concerning the properties of the VEV (|2.10| ) as 



the function of a. It is easy to check that the field (dip) 2 (dip) 2 e 1 ^ has a second order 
resonance with the field e l ( a+2 P' v at a = —(3/2. Similarly, at a = (3/2 it has second 
order resonance with e 1 '-" -2 ^^. Therefore the VEV ( |2.10| ) is expected to have poles at 
a = ±(3/2. As we will argue below the residues at these poles can be expressed through 
the VEV of the primary fields responsible for the resonances. 

To explain this point let us consider a product of two primary fields e iaiV (x) e lOL2V (y) 
in ( |l.lp . The corresponding OPE has the form 

e^(x)e^(y) = [ (r) e^+ n ^ (y) + . . . } , (2.11) 

n= — oo *■ J 

where a = a\ + az, r = \x — y\, and the dots in each term stand for contributions of the 
descendants ( |1 . 2j ) of the field e l< - a+n ^^(y). The coefficient functions C are in principle 
computable within the Conformal Perturbation Theory (CPT) (see also [0). The 
CPT suggests for them the following form 

^n,0 _ ^\n\ r 4a 1 a 2 +4n/3(a 1 +a 2 )+2|n|(l-/3 2 )+2n 2 /3 2 jnfi ^ 2 r ^P 2 ^ ; (2.12) 



where the functions / in ( 2.12 ) admit power series expansions, i.e. 

oo 

/o;°a 2 W = E /rVi,a 2 ) t k . (2.13) 
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The CPT gives the coefficients in (|2.13 ) in terms of certain 2 \n\ + 2 /c-fold Coulomb- 
type integrals. Note that the leading terms /g 1 ' (cti, 01,2) in the series ( |2.13[ ) are expressed 
through the integrals 

1 „ n n n 

j n (a,b, P )=- I n d 2 x k n \x k \ 4a \i-x k \* b n \ Xk - Xp \^ , (2.14) 

k=l fc=l fc<p 



namely 



Jo (ai,a 2 ) = 1 



fo ,0 (a u a 2 ) = j n (a 1 (3,a 2 (3,(3 2 ) for n>0, (2.15) 
f a '°(a 1 ,a 2 ) =jn(-ai(3, -a 2 (3,(3 2 ) for n< . 
The integrals ( |2.14| ) are evaluated explicitly |T3[ , 



j n (a,b,p) =7i n ( 7 (p)) n JJ 7(M X 

fc=l 

71-1 

I J 7(1 + 2a + kp) 7 (1 + 2b + kp) 7 ( - 1 - 2a - 2b - (n - 1 + k)p) 



(2.16) 



A;=0 

Here and below the notation 7 (£) = Y(t)/Y{l — t) is used. Let us quote also the expression 
for the first subleading term in the expansion (|2.12j ) of the function C°'° , 

/°'°(ai, a 2 ) = J(«i/3, « 2 /3, /3 2 ) , (2.17) 

where 

J(a,b,p) = J d 2 xd 2 y \x\ 4a \y\- 4a \l-x\ 4b \l-y\- 4b \x-y\- 4p . (2.18) 

This integral can be expressed through generalized hypergeometric function ^F 2 at unity 
(see also |P0| ). The coefficient functions standing in front of the descendant field in 
( [2. 11|) admit similar CPT expansions. 

There are reasons to believe that the series ( |2.13| ) (and similar CPT series for the coef- 
ficient functions corresponding to the descendant fields in ( |2.11| )) converge for all complex 
t. But independently of the convergence these series can be used to generate asymptotic 
short-distance expansion for the two-point correlation function 

Q aia2 (r) = ( e ia ^(x)e ia ^(y) ) SG , r=\x-y\, (2.19) 

provided the one-point VEV of the exponential fields in the r.h.s. of ( |2.11j ), and also 

the one-point VEV of their descendants ( |1.2j ), are known. For the exponential fields the 
one-point VEV 

G* = {e iatp )sG (2.20) 

closed analytic expression exists ||. According to our discussion above the first nonzero 
contribution due to the descendants comes from the fields (d(p) 2 (B(p) 2 e l ( a + n P)v ; namely 

e« ia2 (r) = £ lc-f a2 (r) ( e *(«+^ ) 5G + (r) ( {dtpf&pfe****) * ) 8G + ■ ■ ■ } , 

(2.21) 
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where a = a.\ + a 2 and the omitted terms contain the descendants of the levels / = / = 4 or 
higher. The coefficient functions C™' 2 a2 (r) admit CPT expansions similar to ( |2.12[ ) , ( |2.13| ) . 
In particular, 

C°f Q2 (r) = r * (l + 0(A 4 - 4 ^ 2 )) ■ (2-22) 

Combining all these expressions one can write down the r — > expansion 



2 

4 



G aia2 (r) = G ai+a2 r 4 ^ 2 jl + J( ai (3,a 2 p,(3 2 ) /x 2 r 4 " 4 ^ - U{pLx + a 2 ) r 

oo 

+ 0(//r 8 - 8 ^ 2 ) } + ^ 1+ a 2+ n/3 jn(«i/3,a 2 /3,/3 2 ) {l + 0( M 2 r 4 - 4 ^)}x 



n=l 

oo 

,2 fl 2 



A* 



n=l 



n 4aia 2 -4n/3(ai+a 2 )+2n(l-/3^)+2n''/3-' 



{l-fO^r 4 " 4 ^)} M ' 

(2.23) 

where TC(a) stands for the ratio 



(dip) 2 (d<p) 2 e ta ^} 



Note that at a = —(3/2 the leading contribution of the field e 1 * 0-1-2 ^^ has the same power 
low in r as the contribution of the descendant (dip) 2 (dif) 2 e' iaip . This is exactly the second 
order resonance we have mentioned above. The contribution comes with the coefficients j 2 
which exhibit the pole at this value of a as is seen from ( |2.16p . The VEV ( |2.10| ) also has a 



resonance pole at this point, and these two pole terms must compensate. This requirement 
leads to the relation 

/ s 2+2? 4 i 

res Q= _^(«) = -^j - (1 + 3 7(" g " 7(0 • (2-25) 

The compensation of the poles at a = —(3/2 results in the logarithmic term r 4a i a a+4 l Q g(r) 
in the short distance expansion of ( p,19|) with a\ + a 2 = —(3/2. The relation ( |2.25| ) will 
be used in the next section to fix the normalization of the VEV ( |2.10|) . 



3. Reflection relations 

In the "reflection property" of the Liouville CFT was used to derive the one-point 
VEV of the exponential fields e iaip in sine-Gordon model. Let us first briefly remind the 
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arguments of 0], and then show how these arguments can be extended to the case of the 
descendant fields (|1.2|). 

The sine-Gordon model ( |1 . 1| ) is closely related to the sinh-Gordon model 

A BhG = j d 2 x ^ -^(d^f + 2/i cosh(M | . (3.1) 

In particular, the one-point VEV (e lotip )sG and (e aip ) s h G in the two models are related 
through the substitution 

b = i/3 , a = ia . (3-2) 



This relation also holds for the one-point VEV of the descendant fields ( |1.2| ). In turn, the 
sinh-Gordon model can be regarded as the Liouville CFT, 

A L = J d 2 x | Y^(^) 2 + fie b< ^ | , (3.3) 

perturbed by the operator e~ b(fi . As is known ||, the correlation functions of the fields e a(p 
in the Liouville theory exhibit the following "reflection property" , 

( e^(x) ...) L = R(a) ( e^- a ^(x) ...) L , (3.4) 

where 

Q = b~ 1 + b (3.5) 

and the coefficient function 

/ ^r(fr 2 ) \ 1+ ^ T(-^ + f)T(-b 2 + 2ab) 
K[a) \r(l-b*)J T(2 + ±-f)r(2 + bi-2ab) [6 - b) 



is essentially the vacuum reflection amplitude of the Liouville CFT. The relation (|3] 
suggest the following "reflection relation" for the one-point VEV of (|3.1| ), 



(e a v) shG = R(a) (e ( - Q - a ^) shG . (3.7) 

Combining this relation with the obvious symmetry property 

(e a *) shG = (e- a *) shG (3.8) 

and certain assumptions about analytic properties of ( e a(p ) s h G one can derive VEV for 
exponential field. 
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Conformal symmetry of the Liouville theory ( |3.3| ) is generated by the energy- 
momentum tensor 

The exponentials e a(p are primary fields with respect to the Virasoro algebra generated by 
flOP- Let us introduce the notation 

L[ n ] £[m] = L— ni L- n2 . . . L— nN L— mi L— m2 . . . L— mK e aip (3.10) 

for the corresponding descendant fields. The symbols [n] and [m] here stand for arbitrary 
strings \—n\, — ti 2 , . . . , — njy], [— mi, -m2, . . . , —ttik]- In ( |3.10|) L n , L n are standard Vira- 



soro generators associated with ( [3.9| ). It is possible to show that the reflection property 
extends to all these descendants, namely 

(L [n] L [m] e at ?(x) ...) L = R(a) (L [n] L [rn] e^- a ^(x) ...) L . (3.11) 



The arguments identical to those in |7J suggest then that the "reflection relation" (|3[ 
generalizes to the descendant fields ( |3.10|) in straightforward way, 

(L [n] L [m] e a *) shG = R(a) ( L [n] L [m] e^~ a ^ ) ahG . (3.12) 



The generalization of the symmetry relation fl3.8p is less straightforward. The relation 
( ^.8| ) is a simple consequence of the symmetry (p — > —99 of the action ( J3.1[ ) . However, while 



the action (|3.1| ) is invariant with respect this transformation, the components ( |3.9| ) of 
the modified energy-momentum tensor, and hence the corresponding Virasoro generators 
L n , L n , are not. In this respect the basis 

(d n V) (d n V) . . . (d nj » (a m V) (9 m V) ■ • • {d mK <p) e atp (3.13) 



in the space of the descendants is more convenient as the fields (|3 .13 ) transform under 



this reflection in an obvious way. The fields ( |3.10|) can be written down as the linear 
combinations of the fields (|3.13|) of the same levels I, I and vice versa @. Finding this 



2 Of course if the Virasoro module with the primary field e aip has a null vector at the level I 
the relation between (3.10) and (3.13) becomes singular. Below we consider generic case of a and 
ignore this subtlety. 
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relation for given levels requires solving a finite algebraic problem, as explained in H. 
Here we will only need the relation ||, 



L_ 2 L_ 2 e a ^ 



\ (^) 2 + (f + a) d\ ) ( - \ (M 2 + (f + a) aV ) e- . (3.14) 



Consider the one-point VEV of the field ( |3.14|) , 

(L_ 2 L_ 2 e^) shG = ^- (l + 2a(Q + 2a)) 2 ( (c^) 2 ^) 2 > sW? 

lb 



(3.15) 



where the relations analogous to ( |2.9|) were used to simplify the r.h.s. Then it follows from 
( P^|) that 



(1 + 2a(Q + 2a)) 2 ( (d^ 2 (d^) 2 ^ ) shG = 



R(a) (1 + 2 (Q - a)(3Q - 2a)) 2 ( (d^) 2 («^) 2 e^~ a ^ ) shG 



(3.16) 



We find that the function 



H(a) 



{d V ) 2 {dipfe^) shG 



shG 



(3.17) 



satisfies the functional equations 



H(a) 



(26 + 3/6 - 2a) (36 + 2/6 - 2a) 



(6 + 2a) (1/6 + 2a) 



H(Q - a) 



(3.18) 



H(a) = H(-a) 



where the second equation follows from the obvious symmetry of (|3.17|) . Note that the 
equation ( |3.18|) remains unchanged if one makes the substitution 



6^ 6" 1 . 



(3.19) 



This is in agreement with well known "duality" symmetry of the sinh-Gordon model ( |3.1|) , 
which in particular implies that all VEV of the fields ( |3.13| ) must be invariant with respect 



to the transformation (|3.19|) . 

Obviously, the equations (|3.18| ) determine the function H(a) only up to a factor F(a) 
which is an even periodic function, 



F(a) = F(-a), F(a)=F(a + Q) . 



(3.20) 
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The solution we are interested in must have the poles at a = ±6/2 corresponding to the 
second order resonances discussed in the previous section. Also, the function H(a) must 
respect the symmetry (|3.19|) . Strictly speaking, this information is not sufficient to fix the 



ambiguity ( |3.20| ). Nevertheless, there is a "minimal" solution which satisfies the above 



requirements, 



H(a) 



8Q 2 



,a b x , a 5 x , a 1 x , a fx 
7 'g ~ 2Q^~ Q ~ 2Q^^Q ~ 2~bQ' 1 ^~Q " 2bQ^ 



(3.21) 

where 7(t) = T(t) /r(f — t) and m is the mass of the sinh-Gordon particle. The residue con- 
dition analogous to (|2.25|) is used to fix the overall normalization of ( |3.21|) . We conjecture 
that this minimal solution gives exact ratio (|3.17|) in the sinh-Gordon model. The VEV 
( p.lO| ) is then obtained by the substitution (pT2|) , which yields ( |1.3| ). In the subsequent 
sections we give some evidence in support of this conjecture. 

4. Comparison to semiclassical results 



The result ( |1.3|) can be checked against certain semiclassical calculations in ( |1 . 1|) . 
Consider the two-point correlation function ( |2.19|) with 

ati=u(3, a 2 = a/f3, (4.1) 

where both a, uj ~ 1, in the limit (3 — > 0. In this limit the functional integral defining ( |2.1 9| ) 
is dominated by the saddle-point configuration (p c i(x) = ^ 4>{t), t = m\x — y\, where (p(t) 
is a solution to the Painleve III equation 

d?(j) + t' 1 d t (f> = - sinh (2(f)) (4.2) 
2 

regular at t > and satisfying the asymptotic conditions 

0(i) = 2a log(8/t) - log ( 7 (± - ct)) + 0(t 2±4 °) as t - , 



, . . 2 sin(7ra) Tr , . 
(j)(t) -> — - K (t) as t -> +00 

TV 



(4.3) 



where Ko(t) is the MacDonald function and again ^y(x) = T(x)/T(l — x). Therefore the 
correlation function under consideration can be written as 



^{x)e^{y)) 



SG 
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= (e 2 <^)~" . (4.4) 



As is known [EH] , this solution to the Painleve III equation admits a double-series expansion 



^ (t) = 7 E {™ + n + 2a(m-n)Y B m , n ~ 



m,n=0 



I \ 2(m+n-l+2cr(m-ri)) 



(4.5) 



where the coefficients B m ^ n satisfy certain recursion relations (see [^TJ for details). Using 
these relations one can derive explicitly few first terms of the expansion (^.5|), 



MM - ^ 



t\ 2 48 B/ l \ /t\4-i2^ 

+ 7 6 (tt+^) 



(l + 2a) 2 V87 (l-2a) 4 
64(1 -2a) o,l N /t\ 4 - 4CT 16 2 ,1 \ /t\4+4 t r 



(l-4(7 2 ) 2 



(4.6) 

This expansion is to be compared with the corresponding limiting case of the expansion 
( |2.23| ). To make this comparison straightforward one can use the the relation 



'a/P+uP+np 



G a /P GojP 



P 2 ^0 



m 



2uj+2n 



(4.7) 



which is obtained from the explicit formula for the VEV (|2.20|) ||, and the following 
limiting expressions for the integrals ( |2.18|) and (|2.16|) 



J{a,ujj3 2 ,l3 2 ' 



(l-4a 2 )< 



n np2n r(2u + n) 



(4.f 



(l + 2a) 2n n\ T(2u) 
Also, assuming (|0|) valid one has for the ratio (|2.24| ) 



TTv 

"w+^)i,.„^ 16j;(1 _ 4 „ r 

Finally, //|^2_, — > ItT^' and with (13)' 0)' ((OD tlie ex P ansion (|2.23| ) takes the form 



(4.9) 



SG 



T(2uj + n 



p 2 ^o 



t \ 4^°" 



1 i 2lj 

7(3 -^). 



t 4 (2(7 + ^1 



0(t 8 ) + £ 



n! r(2w) 



1 + 2(7 



£ \ 2n(l+2(j) 



64 (l-4a 

i + o(t 4 ))+ 



2^2 



(4.10) 



r(-2cj +n) 
" x n\ T(-2w) 



2 7(| + ^) 

1-2(7 



2 n 



£ x 2n(l-2cr) 



i + o(* 4 ) 
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where again t = m\x — y\. It is not difficult to see now that (^.6| ) and (f4.10| ) are in 
exact agreement with (|4.4|) . Thus the semiclassical relation (|4.4| ) actually yields (|4.9| ) 
and therefore supports our main result ( |1-3| ). It is interesting to notice that agreement 
between (|4.10|) and (|4.5| ) further suggests the following explicit expressions for some of the 



coefficients B m>n in (|4.5| ) 



4 n 1 
B ^= n 9 , 2n l 2n {- + a), n=l,2,..., 

n(l — 2a) Zn v 2 ' (4.11' 

-B mj o = -B m+2 ,i = , m = 1, 2, . . . , 



which are not immediately obvious from the recursion relations of |21 . 



5. Comparison to perturbation theory 



It is easy to see that ( |1.3| ) admits power series expansion in /3 2 . For further references, 
let us quote two special cases. First, for a = 



2 4 4 4 

n m m m 



< WW ho = - j-^ = - w + w + Oil) . (54) 

Here we have chosen to express the VEV ( p..6|) through the mass m = 2Msin(7r^/2) of the 
lightest sine-Gordon breather. Second, if a = u(3 where u is a constant, 

These expansions can be compared with the results obtained directly from ( |1 . 1| ) by means 
of ordinary Feynmann perturbation theory. 

In developing the perturbation theory it is convenient to start from the action ( |1 . 1] ) 
in its "bare" form 



ASG = J d 2 X ^{ \{d^f - ^ [00B(^)] fl } = 



(5.3) 



where the symbol [■■■]b signifies that we are dealing with the bare fields (as opposed to 
the renormalized fields defined in Sect.2), and it is assumed that the action ( |5.3|) is supple- 
mented with some cutoff procedure, with the cutoff momentum A. The scale dimensions of 
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the bare fields coincide with their naive values and hence the bare mass parameter mo has 
the dimension of mass. The relation between mo and the physical mass m of the lightest 



breather particle of ( |1 . 1| ) can be found perturbatively, order by order in /3 2 . For instance 



with the account of the leading mass correction diagram in Fig.l, we have 

ml=m 2 + m 2 l3 2 L + 0(/3 4 ) , (5.4) 

where 

r 1 f d 2 k ,A 2 ^ 

n J k -\- m z m z 
and C is a constant whose exact value depends on the implementation of the cutoff pro- 
cedure. In fact, to all orders in f3 2 this relation has the form ||10|| il 

2 



m = m 



2 e f3L h(f3 2 ); h([3 2 ) = l + —[3 4 + 0((3 G ). (5.6) 



6 



Similarly, the relation between the bare exponential fields [e 10 ^] B and corresponding renor- 
malized fields can be written as 



a 2 

(~AV^) [e^] fl , (5.7) 



where 7 is Euler's constant and the normalization of the field e %OLip is fixed by the short- 
distance asymptotic condition 

(e^(x) e-^(y)) SG ^\x-y\- 4a2 as \x - y\ ^ . (5.8) 

Now we are prepared to make some perturbative calculations of the one-point VEV. 
As the first example let us consider the field (dp) 2 (dp) 2 . According to our discussion 
in Sect. 2 this renormalized field differs from the corresponding bare field by appropriate 
counterterms 

(dip) 2 (8ip) 2 = [(dip) 2 (dip) 2 ] B + m 4 (A + A 1 [cos(2(3p)} B )+m 2 A 2 dd[co S ((3p)} B , (5.9) 



where Aq,A\,A2 are constants which can depend on (3. In ( |5.9| ) the counterterms of the 
type A 2 [c^d^] and A 4 which are needed to absorb the quadratic divergences in the 
matrix elements of [(d(p) 2 (5(p) 2 ~\ B are not written down. In the following calculations of 
these matrix elements we will systematically subtract all quadratic divergences; with this 



Exact form of the function h(/3 ) can be found in [ 10 1 . 
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convention the "quadratic" counterterms can be ignored altogether. The counterterms 
explicitly shown in ( |5.9[ ) are to compensate for remaining logarithmic divergences. It is 
possible to see that this compensation can not be achieved with the coefficient A2 being just 
constant; instead one has to set A2 = AL + B, where L is the logarithm ( |5.5| ). The reason 
for this subtlety lays in the fact that the field {dip) 2 {dip) 2 always has a first order resonance 
with the field dd cos{f3ip) which results in the logarithmic scaling of all its matrix elements 
which receive contributions from the above total derivative field. Fortunately, here we are 
interested only in the one-point VEV 

( {dip) 2 {dip) 2 ) SG = { [{d^ 2 {dip) 2 ] B ) SG + A m^ L h 2 {(3 2 ) + A im ^e- 2 ^ L h 2 {(3 2 ) g{2$) , 

(5.10) 

which gets no contribution from the last counterterm in ( |5.9| ) and hence is not sensitive 
to the above subtlety. In writing ( p.lO| ) we have used ( |5.6| ) to express mo through the 
physical mass and also used the notation 

{[e^] B ) SG = e- a2L g{a) . (5.11) 

Explicit expression for g{a) can be found in ||; here we will only use the fact that for 
fixed oj 

9 {u;f3) = l + 0{f3 6 ) . (5.12) 

The first term in Q5.10] ) can be calculated directly using Feynmann diagrams for ( |5.3| ). To 
the leading order in f3 2 one obtains 

( [{dip) 2 {dip) 2 ] b )sg = \ m*L 2 + 0{(3 2 ) . (5.13) 

(Let us remind that we subtract the quadratic divergences). In order to compensate for 
this 1? divergence the counterterm coefficients in ( |5.1U|) have to be chosen as follows 

Ao = -^ + 0(p, A 1 = -^L + 0(^) (5.14) 

and we obtain 

({d^ 2 {dp>) 2 )sG = -j^ + 0{^) . (5.15) 

The calculation can be easily extended to the next order in /3 2 . The next perturbative 
contribution to the VEV ( |5.13| ) comes from the diagram in Fig. 2. It has the form 

m 4 (3 2 (aiL 2 + a 2 L + a 3 ) , (5.16) 
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where a-\_,a,2, a 3 are numerical coefficients. It is not difficult to check that in order to find 
the next term in ( |5.1|) one only needs to know the coefficient a\ in front of the leading 
logarithmic term in (|5.16| ). This coefficient is evaluated directly from the diagram, a± = 1. 
Then compensation of this term requires the following terms in the f3 2 expansion ( |5.14| ) 

A ° = -w + w +om - Ai = -w + w +om - (5 ' 17) 

The finite terms remaining in (|5.10| ) after the cancelation of the divergences yield exactly 
(0)- 

Next, let us apply the perturbation theory to more general VEV ( [2.10|) with a ^ 0. 
Again, the renormalized field is a combination of corresponding bare field and suitable 
counterterms, 

(A 2 e 2 7+ C /4 p 2 (0^)2(^)2^= [(a^) 2 (^) 2 e^] B + 

il dd(B+ [e^ a+ ^] B +B_ [e^ a -^] B )+ (5.18) 



4 



The constants A, B have to be determined from the requirement that the renormalized 
field {dip) 2 {dip) 2 e iaip has definite scale dimension 4 + 2a 2 . It is convenient to divide ( |5.18| ) 
by the VEV ( |5.11|) , and trade the parameter mo in favor of m, 



(dip) 2 (dip) 2 _ [(dip) 2 (5ip) 2 e^] B ^ 



(e^)sG ([e ia v] B )sG 
m 2 h((3 2 ) dd(B + e- 2 ^ L 9 ^ a + 0) [ e '(«+«*] N + B. e 2 ^ L 9{a ~ 0) [e^-^] N ) + 

m 4 h 2 ((3 2 ) (A + e~W 2+4 ^ L ^ ± W) \ e ^+ 2 ^} w + A e 2 ^ L + 
V 9\p) 



(5.19) 



A _ e -(2/3 2 -4a/3)L 9{<* ~ ri( a -2/3) v l 

g(a) iN 
where we used the notation 

[e iaip ] N = e^/(e^)sG=[e^] B /([e^] B } SG . (5.20) 

One can calculate perturbatively matrix elements of ( |5.19j ) , adjusting the coefficients 
order by order in (3 to ensure the cancelation of all L-dependent terms. Note that ( |5.20j ) 

contains no divergences and so all the L dependence of the counterterm part in ( |5.19|) is 
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shown explicitly. It turns out that contrary to the case a = studying just the VEV of 
( |5.19|) is not enough to determine the coefficients A+, Aq, A-. We have considered the 
matrix elements of (|5.19| ) between the vacuum and one- and two-particle states (involving 



the lightest breather) along with the VEV. The calculations are straightforward but rather 
bulky and we do not present them here. In the case a = u(3, 10 ~ 1 and in the leading 
order in f3 2 the cancelation of L-dependent terms requires the following choice of the 
coefficients, 

B+ = W [ 2u(l+u) 2 + °^' B ~ = 2u(l-u) 2 + ' 

and 

1 1 n(—\ A - -i - o(— \ 

+ " (3 4 16(l + o;)(l + 2a;) ' ~ ^ 16(1 - w)(l - 2w) + ^ f3 2 > ' 

A ° = 'W 8(l + w)(l-u;) + • 

(5.21) 

With Q5.21 ) the result for the VEV of this field identical to ( |5.2|) immediately follows from 

(EH)- 



6. Exact results from XYZ model 



As is well known [R2| , the sine-Gordon QFT (|1 . 1|) can be obtained by taking an 
appropriate scaling limit of the XYZ spin chain described by the Hamiltonian 

1 N 

Kxyz = - Ye ( J * < + J v < ° v s+i + J * < <+i - J ) > (6- 1 ) 

s=l 

with J x > J y > \J Z \. In ( |6.1| ) we have introduced an auxiliary parameter e which is inter- 
preted as a lattice spacing. It is convenient to use the Baxter's elliptic parameterization 
P^j of the coefficients J in (|Q|), 

J 



tt V^(o)^i(/3 2 ) 4 (o)0 4 (/3 2 ) 
i-p 2 fe 4 ((3 2 )e[(o) d x (o)^(o) 



n V^(o)^i(/3 2 ) e 4 (o)e 4 (p 2 

i - p 2 { sap 2 ) m 



2\ 



j_ 1 - P 2 / 6[((3 2 ) | 6' 4 (P 2 



n V<9i(/5 2 ) W 2 ) 
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where 

oo 

d 1 (v) = 2p* sm(nv) JJ (l - p 2n ) (l - e 2niv p 2n ) (l - e" 2 ^> 2n ) , 

71=1 

OO 

9 4 (v) = Y[ {l-p 2n ) (1 - e 2wiv p 2n ~ x ) (1 - e ~ 2 ™ v p 2 ™- 1 ) 

n=l 

and the prime in (|6.2| ) denotes a derivative. The scaling limit of ( |6.1| ) is achieved by sending 

iV^oo, e^O, p^O (6.3) 

with the combinations 

R = Ne, M = ^p (1+ ^/ 4 (6.4) 

kept fixed. According to Refs. p3| , |2"2"| ] in this limit the energy spectrum of ( |6.1[ ) is described 
by the QFT ( p. . 1| ) , the parameter M coinciding with the sine-Gordon soliton mass. 

In fact, the QFT ( |1.1| ) itself controls only the leading p —>■ singularities in the spec- 
trum of (|6.1|) . Using exact XYZ ground state energy |2cj] one can easily extract subleading 
singular terms in this quantity. Being expressed through the scaling parameters M and R, 
the singular at p — ► part of the bulk ground state energy reads 

{E X Yz) sl n g = -^f- tan«/2) j 1 + (^ff + 0(e 4 ) J (i?M»l). (6.5) 

Whereas the leading term here is exact sine-Gordon vacuum energy, the higher-order in e 
terms must be attributed to the irrelevant operators which differ the exact XYZ Hamilto- 
nian (|6.1|) from the Hamiltonian Hsg of the sine-Gordon QFT ( |1.1|) . As follows from the 
analysis in ||25|| , for f3 2 < 2/3 the leading in e correction comes from the terms 

U XYZ = const + U SG -- ^ — (a+(^) 2 (^) 2 + A_((^) 4 + (^) 4 ))+... , (6.6) 

where A_|_ and A_ are numerical coefficients whose exact values are found in |25|] and the 
dots stand for the irrelevant operators of higher dimensions. The corrections in ( |6.5| ) can 
be expressed through the expectation values of the correction terms in ( |6.6| ) over the sine- 
Gordon vacuum. Obviously, it is the VEV of the operator [dip) 2 (dip) 2 which is responsible 
for the e 2 term in ( |6.5| ) (the VEV of (dip) 41 and (d(p) 4 vanish), i.e. 

M 4 

A + ( (dp) 2 (dp) 2 ) SG = — tanK/2) . (6.7) 



Using the result of [J25 



cotK/2) 

K = 7, > ( 6 - 8 ) 

one arrives precisely at the Eq.( |1.6| ) which therefore agrees with exact results of the lattice 
theory. 



19 



7. Application: Two-point correlation function in scaling Lie- Yang model 



In Ref. a two-point correlation function in so called Scaling Lee- Yang Model (SLYM) 
was studied. In particular, a combination of the operator product expansions and confor- 
mal perturbation theory was used there to develop a short-distance expansion similar to 
( |2.23| ). In this section we will use our result to extend this expansion further thus 

obtaining more accurate estimate for the two-point correlation function at all distances. 

The SLYM is one of the simplest of the perturbed CFT (|1.7|) , namely 



Aslym = M 2 /s + ih J d 2 x $(x) , (7.1) 



where 

= $i, 3 (x) , A$ = -- . 

5 

As is known (see e.g. [|15|j ) the QFT ( f7.1|) is massive; it has one sort of massive particles 
whose mass m is related to the parameter h in ( |7.1|) as 

2*5t (r(2/3)r(5/6))^ 12 12 . . 

h = T~r — w mwr// m ~ = °- 0970485 ---m- . 7.2 
16tt5 r(3/5)r(4/5) 

We will use the notations 

9(x) = TZ(x)/4 = ihn (1 - A*) $(x) (7.3) 

for the trace of the energy- momentum tensor associated with ( |7.1| ). 
Consider the two-point correlation function 

G(t) = (e(x)e(O)) , r=\x\. (7.4) 

According to || this correlation function admits the following short-distance expansion 

G(r) = - h 2 n 2 (1 - A,,) 2 C^(r) + ihn (1 - A,,) C* $ (r) ( G}- 

/iV(l-A„) 2 Cllir) (Tf)+0{rf) , ' ' 



where the notation 

Tf = L_ 2 L_ 2 I (7.6) 
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is used. The coefficient functions C in (73) admit power series expansions in h, the first 
few terms being known expiicitly Q @ 

C 1 (r) -rUl + ^- r4(1/5)F(3/5) hr% +0(r*) \ 
6w(r)_r l 1 + 1960 rs(4/5) 

c * (r) _ i r 2 (l/5)r(2/5) g r 2^ I*(l/5) I«(2/5) f f , 

Gw(r) " * IOtt T(4/5) r + 5^ 9 r3(3/5) 1^(4/5) hr + ^ > J ' 



2 4 



C*^(r) = ^{l + 0(^)} 



(7.7) 

With the known exact VEV of the field 0, 

(6) = - i ^=m 2 , (7.8) 

Eqs. (|7.5|) , (|7.7[) effectively gives the short-distance expansion of the correlation function 
( [Tjj] ) up to the terms ~ r 2 ^ [EI. Now, using ( |1.8| ) we can derive the VEV 

(TT) = -^m 4 . (7.9) 

24 

This additional peace of data allows one to compute explicitly the term ~ r~ in (|7.5|). 
The next term ~ r^, which would come from the h 2 term in Cj^,, is still not available in 
an analytic form. 

The correlation function (|7.4|) admits also the large-distance expansion in terms of 
exact form factors 0. Two leading terms, corresponding to zero- and one-particle contri- 
butions, are known in analytic form, 

G(r) = ^m*{l-^ ZK (mr) + ... }, (7.10) 



wV3tt r [— dt t 

exp \ - / — - } = 0.8155740 



where 

27 I ,/ 7r sm(i) 

and -RTo(t) is the MacDonald function. Further terms in this expansion of G(r) can be 
obtained by numerical integration of its spectral representation including the contributions 
of two or more particles in the intermediate state || . The expansion is known to converge 



4 Notice the analytic expression for the first subleading term in the expansion of C*<i>, which 
was given numerically in [El. 
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very fast. With the inclusion of up to four-particle contributions this expansion gives a 
precision better then 10 _2 % for mr > 10 -2 . The short-distance expansion ( |7.5| ) (with 
( [7.7|) , ( |7.8|) and (|7.9|) ) is compared with this data in Table 1. The combined data from 
these two expansions apparently have relative precision 10 _5 % or better for all values of 
r. 

Finally let us note that since exact form factors of the sine-Gordon model are known 
]2^,^|, similar numerical analysis can be performed for the general sine-Gordon correlation 
function pJ9|) . 
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Long-distance expansion 


Short-distance expansion 


mr 


0-1-2-3-4 particles 


Without (TT) 


With ( TT ) 


0.001 


0.01964182 


0.01947405 


0.01947405 


0.002 


0.02553233 


0.02547150 


0.02547150 


0.005 


0.03616744 


0.03615515 


0.03615515 


0.010 


0.04689571 


0.04689286 


0.04689286 


0.020 


0.06045855 


0.06045806 


0.06045806 


0.040 


0.07730747 


0.07730741 


0.07730741 


0.060 


0.08881026 


0.08881025 


0.08881025 


0.080 


0.09771443 


0.09771443 


0.09771443 


0.100 


0.10502922 


0.10502922 


0.10502922 


0.120 


0.11125352 


0.11125352 


0.11125352 


0.140 


0.11667524 


0.11667523 


0.11667525 


0.160 


0.12147752 


0.12147749 


0.12147752 


0.180 


0.12578494 


0.12578489 


0.12578495 


0.200 


0.12968660 


0.12968651 


0.12968661 


0.220 


0.13324863 


0.13324850 


0.13324866 


0.240 


0.13652170 


0.13652150 


0.13652174 


0.260 


0.13954552 


0.13954523 


0.13954559 


0.280 


0.14235191 


0.14235151 


0.14235201 


0.300 


0.14496678 


0.14496622 


0.14496693 


0.400 


0.15581014 


0.15580799 


0.15581079 


0.500 


0.16401350 


0.16400745 


0.16401562 


0.600 


0.17045818 


0.17044413 


0.17046372 


0.700 


0.17564907 


0.17562051 


0.17566157 


0.800 


0.17990480 


0.17985216 


0.17993010 


0.900 


0.18344031 


0.18335022 


0.18348740 


1.000 


0.18640771 


0.18626233 


0.18648980 


1.200 


0.19105758 


0.19072646 


0.19127219 


1.400 


0.19446608 


0.19380550 


0.19494924 


1.600 


0.19700864 


0.19581239 


0.19798354 


1.800 


0.19892990 


0.19691763 


0.20073902 


2.000 


0.20039614 


0.19720202 


0.20353862 


2.500 


0.20275757 


0.19434989 


0.21284363 


3.000 


0.20402333 


0.18570096 


0.23007149 



Table 1. Comparison of short and long-distance expansions for the two-point 
correlation function ( [7.4j) . The first column gives the results of long-distance 
expansion which includes contributions of up to four-particle states (the four- 
particle contribution which we include here represents the improvement over the 
data in []3| ). The data in the second and the third columns correspond to the short 
distance expansion ( [7.5| ) without the TT term and with this term, respectively. 
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Figures 





Fig. 1 



Fig. 2 



Fig.l. The leading mass correction diagram which gives ( |5.4j ) . 

Fig. 2. Diagram contributing to the VEV ( |5.13|) in the order 1 . The strokes over the 
propagators stand for the derivatives. 
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